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Based on the first-principles electronic structure calculations, we predict that a chromium oxide 
K2CrgOi6 of hollandite type should be a half-metallic ferromagnet where the Fermi level crosses 
only the majority-spin band, whereas the minority-spin band has a semiconducting gap. We show 
that the double-exchange mechanism is responsible for the observed saturated ferromagnet ism. We 
discuss possible scenarios of the metal-insulator transition observed at low temperature and we 
argue that the formation of the incommensurate, long-wavelength density wave of spinless fermions 
caused by the Fermi-surface nesting may be the origin of the opening of the charge gap. 
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I. INTRODUCTION 



puzzle. 



Half-metallic ferromagnets 1 offer a unique opportu- 
nity for studying the electronic states of strongly cor- 
related electron systems. Here, only the majority- 
spin electrons form the Fermi surface with a gapped 
minority-spin band^ and can couple with excitations of 
the spin (and in some cases orbital) degrees of freedom 
of the system. The situation therefore should attract 
much interest, in particular, when the relevant electrons 
are strongly correlated, leading the system to double- 
exchange ferromagnetism 3 and metal-insulator transition 
(MIT). In this paper, we will show that a chromium oxide 
K2Cr80i6 with the hollandite- type crystal structure^ 
belongs to this class of materials and can provide a good 
opportunity for further development of the physics of 
strong electron correlations. 

The crystal structure of K^CrgOie (see Fig. 1) be- 
longs to a group of hollandite-type phases where one- 
dimensional (ID) double strings of edge-shared Cr06 oc- 
tahedra forms a CrgOi6 framework of a tunnel structure, 
wherein K ions reside. 6 Cr ions are in the mixed-valent 
state of Cr 4+ (d 2 ) : Cr 3+ (d 3 ) = 3 : 1, and hence with 
2.25 electrons per Cr ion. It has recently been reported 7 
that the phase transition from the Pauli-paramagnetic 
metal to ferromagnetic metal occurs at T c ~ 180 K by 
lowering temperatures, where the ferromagnetic state has 
a full spin-polarization of 18 /ie per formula unit (f.u.) 
at low temperatures, which is a realization of saturated 
ferromagnetism. 

In addition to this phase transition, it has also been 
reported 7 - that another phase transition occurs from the 
ferromagnetic metal to ferromagnetic insulator at Tmi — 
80 K, suggesting that the charge gap opens below Tmi- 
Surprisingly, the spin polarization is hardly affected by 
this MIT and no structural distortions associated with 
this MIT have been observed so far* 7 - The mechanism 
of the MIT of this material has therefore been a great 




FIG. 1: (Color online) Schematic representation of (a) 
the unit cell of the body-centered tetragonal lattice (solid 
lines) and (b) Brillouin zone of K^CrsOie. In (a), 
the primitive unit cell is also shown in the thin dotted 
lines. In (b), the symbols represent r(0,0,0), M(27r/a, 0, 0), 
X(7r/a, 7r/a, 0), P(jr/a, ir/a, 7r/c), Ki(0, 0, 7r(l/c + c/a 2 )), and 
K2 (2ir/a, 0, 7r(l/c — c/a 2 )), where Ki and K2 are equivalent. 



In this paper, we perform the first-principles electronic 
structure calculations based on the generalized gradient 
approximation (GGA) in the density-functional theory 
(DFT) and we predict that the materials ^CrgOie (A 
= K and Rb) belong to a new class of half-metallic fer- 
romagnets, i.e., the majority-spin electrons are metallic, 
whereas the minority-spin electrons are semiconducting 
with a band gap. We also show from the GGA and 
GGA+Z7 calculations that the double-exchange mecha- 
nism is responsible for the observed saturated ferromag- 
netism. We then discuss possible mechanisms of the MIT 
of this material and argue that the formation of an in- 
commensurate, long- wavelength spin and charge density 
wave (DW) due to Fermi-surface nesting may be the ori- 
gin of the MIT of this material. The opening of the gap 
in the majority-spin band should, however, be detrimen- 
tal to the double-exchange ferromagnetism, for which we 
consider possible reconciliations. We hope that our re- 
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suits will encourage further experimental studies of this 
intriguing material. 

This paper is organized as follows: In Sec. II, we 
show our method of calculations. In Sec. Ill, we present 
our results of calculations and discuss the origin of fer- 
romagnetism and possible mechanisms of the MIT of 
K^CrgOie. Summary of our work and prospects for fu- 
ture studies are given in Sec. IV. 



II. METHOD OF CALCULATION 

The electronic structure calculation in GGA is car- 
ried out by employing the computer code WIEN2k, 8 
which is based on the full-potential linearized augmented- 
plane-wave (FLAPW) method. We use the exchange- 
correlation potential of Ref. 9 . The spin polarization is 
allowed. The spin-orbit interaction is not taken into ac- 
count. We use 1,221 k points in the irreducible part of 
the Brillouin zone in the self-consistent calculations. We 
use the plane-wave cutoff of K max = 4.24 Bohr" 1 . The 
GGA+Z7 calculation 10 is also made to see the effects of 
on-site electron correlation U on the band structure. We 
assume the experimental crystal structure of K^CrgOie 
observed at room temperature with the lattice constants 
of a = 9.7627 and c = 2.9347 A. 6 The Bravais lattice 
is body-centered tetragonal and the primitive unit cell 
(u.c.) contains four Cr ions, one K ion, and eight O 
ions, i.e., KC^Og. The local structure around Cr ions 
is shown in Fig. 2. We use the code XCrySDenii for 
graphical purposes. 




FIG. 2: (Color online) Schematic representations of the local 
structure of K^CrsOie, together with the three tig orbitals 
d xy , d yz , and d zx of Cr ions in the xyz-coordinate system. 
There are two inequivalent O sites, 0(1) and 0(2). All the 
Cr sites are equivalent. 



III. RESULTS AND DISCUSSION 
A. Ground state 

The total ground-state energies obtained in GGA are 

— 10816.1918 Ryd/u.c. in the ferromagnetic state and 

— 10815.9597 Ryd/u.c. in the paramagnetic state, result- 
ing in the energy gain of 3.16 eV/u.c. due to stabiliza- 
tion by the spin polarization, which is realized already 
in the GGA level. In the ground state, we have the 
full spin polarization with the magnetic moment of 9.000 
/iB/u.c., which consists of the contributions from Cr ions, 
2.043 /i B /atom, 0(1) ions, 0.014 /i B /atom, and 0(2) 
ions, —0.082 /iB/atom, in the atomic spheres, and 1.098 
//b/u.c., in the interstitial region. Note that the 0(2) 
ions inside the double strings of Cr ions have a negative 
spin polarization, whereas the 0(1) ions connecting the 
double strings have a positive spin polarization. 



B. Density of states 

The calculated density of states (DOS) is shown in 
Fig. 3 in a wide energy range covering over the O 2p 
and Cr 3d bands. We find three separate peaks in both 
the majority and minority spin bands: in the majority 
(minority) spin band, the O 2p weight is located mainly 
at -7.5 < e < -1.8 (-7.1 < e < -1.8) eV, Cr 3d weight 
with the t 2g symmetry at -1.7 < e < 0.7 (0.7 < e < 2.6) 
eV, and Cr 3d weight with the e g symmetry at 1.7 < e < 
4.0 (2.9 < e < 5.1) eV. The hybridization between the 
O 2p and Cr 3d bands is significantly large. The Fermi 
level is located at a deep valley of the t^g majority-spin 
band while it is located in the energy gap between the 
O 2p and Cr 3d t^g bands in the minority-spin band. 
Thus, the half metallicity of this material is evident in 
the calculated DOS. 

The calculated orbital-decomposed partial DOS, p a (e) 
(a = xy^yz^zx\ in the Cr 3d t2 g region are shown in 
Fig. 4, where the two components are exactly degenerate, 
Pyz(z) = Pzx{z), for both paramagnetic and ferromag- 
netic states. The three t<i g orbitals are almost equally 
occupied by electrons in the paramagnetic state. In the 
ferromagnetic state, however, the d xy orbitals is almost 
fully occupied by electrons and therefore holes are only 
in the d yz and d zx orbitals. Also, the d xy component has 
a rather high peak-like structure at ~0.7 eV below the 
Fermi level, indicating an essentially localized character 
of the d xy electrons. The d yz and d zx components, on 
the other hand, have a rather wide band spreading over 
—0.2 < e < 0.7 eV around the Fermi level, indicating an 
itinerant character of the d yz and d zx electrons. Here, 
the admixture of the 2p z state of 0(2) is significantly 
large: it occurs between the 2p z orbitals of 0(2) and 
the d yz and d zx orbitals of its neighboring Cr ions (see 
Fig. 2). These results suggest that the double-exchange 
mechanism 3 should work for the occurrence of ferromag- 
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FIG. 3: (Color online) Calculated result for the DOS (per 
spin per formula unit (f.u.)) of K^CrsOie in a wide energy 
range for the ferromagnetic state. The Fermi level is set to 
be the origin of energy. 



netism, as will be discussed further below. 
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FIG. 4: (Color online) Calculated orbital-decomposed DOS 
(per spin per formula unit (f.u.)) of K^CrsOie near the Fermi 
level for the (a) paramagnetic state and (b) ferromagnetic 
state. The Fermi level is indicated by the vertical line. Con- 
tributions from the d yz and d zx orbitals (thin dashed and 
dotted lines) are exactly degenerate. 



C. Band dispersion 

The situation may be clarified further if one observes 
the calculated band dispersion near the Fermi level. The 
result is shown in Fig. 5. We find that a rather disper- 
sionless narrow band of predominantly d xy character is 
located at ~0.7 eV below the Fermi level, extending over 
a large region of the Brillouin zone. On the other hand, 
the dispersive t<i g bands of predominantly d yz and d zx 
character with strong admixture of the 2p z state of 0(2) 
are located at —0.2 < e < 0.7 eV and cross the Fermi 
level. Thus, we have the dualistic situation where the 
essentially localized d xy electrons at ~0.7 eV below the 
Fermi level interact with the itinerant d yz and d zx elec- 
trons of the bandwidth comparable with the intraatomic 
exchange energy of ~1 eV, whereby the Hund's rule cou- 
pling gives rise to the ferromagnetic spin polarization via 
the double-exchange mechanism. 3 

To support this further, we make the GGA+Z7 calcu- 
lation for the present material (of which the results are 
not shown here). We find that, as U increases, the d xy 
band shifts further away from the Fermi level, leaving 
essentially no weight above the Fermi level, whereas the 
d yz and d zx bands with strong admixture of the 0(2) 2p z 



states are much less affected by the presence of U. These 
results are consistent with what is expected in the double- 
exchange mechanism of ferromagnetism. We point out 
that the situation is very similar to the case of Cr02, 
which has so far been discussed in detai h 12 i 13 i 14 

Another aspect noted in Fig. 5 is that we have the 
semi-metallic band structure in the sense that the num- 
ber of electrons is equal to the number of holes near the 
Fermi level (also see Fig. 6 below); there are 12 bands 
for the t2 g orbitals in the unit cell, which are fully spin- 
polarized and are occupied by 9 "spinless fermions", i.e., 
by 9 up-spin (or majority-spin) electrons. The integer 
filling of spinless fermions can give rise to either semi- 
conducting or semimetallic band structure, of which the 
latter is realized in the present system. This result sug- 
gests that the MIT observed in K^CrgOie can possibly be 
a semimetal-to-semiconductor transition, where the band 
gap opens between the lower 9 bands and upper 3 bands 
due to any unknown reasons. To realize this, however, 
we need the symmetry reduction (due to lattice distor- 
tion) because the band degeneracy between the third and 
fourth bands (counted from the top of the 12 t^g bands) 
at the P point of the Brillouin zone cannot be lifted with- 
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out breaking the 4-fold rotational symmetry around the 
c axis of the tetragonal lattice. Without the lattice dis- 
tortion, we could have a zero-gap semiconductor, but the 
experimental data on the temperature dependence of the 
electrical resistivity 7 suggest that this possibility should 
be ruled out. Also, because no indications of this type 
of lattice distortions at the MIT have been observed so 
far j 7 - the semimetal-to-semiconductor transition seems to 
be unrealistic in K 2 Cr 8 0i6. 
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FIG. 5: (Color online) Calculated majority-spin band disper- 
sion of K2CrsOi6 near the Fermi level. There are 12 t^ g bands 
of the Cr 3d orbitals, 3 of which cross the Fermi level. The 
line width is in proportion to the weight of the d xy component 
of Cr in (a) and to that of the 2p z component of 0(2) in (b). 



D. Fermi surface 

The calculated Fermi surfaces of K^CrgOie in the fer- 
romagnetic state are shown in Fig. 6. There are 12 t<i g 
bands, 3 of which cross the Fermi level and form the 
semimetallic Fermi surfaces; i.e., the second and third 
bands (counted from the top) form the electron Fermi 
surfaces (see Figs. 6 (b) and (c)) and the fourth band 
forms the hole Fermi surface (see Fig. 6(a)). The wave 
functions at the Fermi surfaces have predominantly d yz 
and d zx character with large admixture of the 0(2) 2p z 
states as shown above. 

We find in Fig. 6(a) that there is a pair of the ID- 
like parallel Fermi surfaces, which are seen to have a 
very good nesting feature. The nesting vector is aligned 
roughly along the T-Ki direction and has the value q* ~ 



(0,0,0.147)2tt/c or (0, 0, 0.853)2tt/c (a deviation in the 
(q x ,Qy) component will be discussed below). Thus, the 
Fermi-surface instability corresponding to the wavenum- 
ber </*, leading to formation of the incommensurate, long- 
wavelength (with a period of ~7c in the real space) DW, 
may be relevant with the opening of the charge gap in 
the present material. Note that the spin and charge DWs 
occur simultaneously with the same wavenumber qr* since 
we have only the up-spin electrons. 




FIG. 6: (Color online) Calculated Fermi surfaces of K^CrsOie 
in the ferromagnetic state. The 61st to 63rd bands counted 
from the lowest are shown in (a) to (c), respectively. 



E. Generalized susceptibility 

To confirm the nesting features more precisely, we cal- 
culate the generalized susceptibility defined as Xo(q) — 
Ek (f( E k) - f(E k+q ))/(E k+q - E k )^ where E k is the 
band dispersion and / is the Fermi distribution function 
at T = K. The calculated result is shown in Fig. 7, 
where the contribution from only the 61st band is given. 
Contributions from other bands including interband con- 
tributions are rather monotone functions of k over the en- 
tire Brillouin zone. Thus, the peak structure coming from 
the 61st band remains as a maximum even in the total 
susceptibility estimated in the constant matrix-element 
approximation. 

To be more precise, the peak structure at q* = 
0.2957r/c and 1.7057r/c seen in Fig. 7 remains strong, ir- 
respective of the value of (q x •>%)•> although there is a 
small variation in the (q x ,q y ) plane. The true maximum 
appears at ~ (fr/a, 7r/a, g*), or around (ir/a, ir/a, g*), 
splitting and deviating slightly from (7r/a, 7r/a, q*) in the 
(QxiQy) plane. Thus, if we include the effects of elec- 
tron correlations, the susceptibility can diverge at this 
momentum qr*, resulting in the formation of the incom- 
mensurate, long-wavelength charge and spin DW, which 
we hope will be checked by experiment in near future. 
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FIG. 7: (Color online) Generalized susceptibility Xo(q) for 
the noninteracting band structure of K^CrgOie- Contribution 
from only the 61st band is shown. 



F. Metal-insulator transition 

Now, let us discuss possible mechanisms of the MIT 
observed in K^CrgOie at T c ~ 90 K. There are two types 
of scenarios: (i) correlation scenario where the electron 
correlations play an essential role and (ii) band scenario 
where the band structure solely determines whether the 
system is metallic or insulating. As a correlation sce- 
nario, we may have a possibility of the formation of the 
incommensurate, long- wavelength charge and spin DW 
due to the Fermi-surface nesting as discussed above. This 
may give rise to the opening of the charge gap, result- 
ing in the MIT observed in the present material. As 
another correlation scenario, we may have the charge or- 
dering (CO) if the intersite Coulombic repulsions between 
charge carriers are sufficiently strong. In the present sys- 
tem, there may occur the localization of 3 spinless t<i g 
holes per 4 Cr site. We find that the possible spatial 
CO patterns are equivalent to those of K2V80i6 j 15ll6llT 
which results in the doubling of the unit cell along the c 
axis. As a band scenario, we may have the semimetal-to- 
semiconductor transition in the present system with the 
integer filling of spinless fermions as discussed in Sec. Ill 
C. Here, no spatial localization of carriers occurs, but 
the band gap may open if the 4-fold rotational symmetry 
around the c axis is broken due to lattice distortions and 
the distortion is sufficiently large. 

Although the lattice distortion is essential for the MIT 
in any of the above cases and its identification can spec- 
ify which scenario is realized, no structural distortions 
associated with the metal-insulator transition have been 
observed so far. 7 This seems to suggest that the simple 
distortions such as the doubling of the unit cell along 
the c axis and the breaking of the 4-fold rotational sym- 
metry around the c axis may be ruled out. We there- 
fore argue that the formation of the incommensurate, 
long-wavelength DW, the observation of which is some- 
times not very easy, may be relevant with the MIT in 
the present material. We hope that the anomaly at 



estimated above will be sought for carefully in future ex- 
periment. 

Another experimental fact to be noted is that the spin 
polarization is hardly affected by the MIT in this ma- 
terial. This is interesting because the double-exchange 
mechanism for ferromagnetism (of which the energy gain 
occurs due to the first-order process of the hopping) 
should cease to work if the charge gap opens clearly in 
the majority-spin band and coherent motion of carriers 
vanishes. The double-exchange ferromagnetism may be 
killed in such cases. However, if the charge gap is small 
enough, the second-order processes of the hopping of con- 
duction electrons can give rise to the ferromagnetic in- 
teraction between localized spins as well and thus the 
ferromagnetism can be maintained. The vanishing first- 
order process seems to affect very little on the already 
fully spin-polarized system. Experimentally, 7 the elec- 
trical resistivity measurement shows that the tempera- 
ture dependence is of a three-dimensional variable-range- 
hopping type rather than a thermal activation type and 
thus it is not clear whether the well-defined charge gap is 
present in this material. Also, the charge gap expected 
in the formation of the incommensurate, long- wavelength 
DW state corresponding to q* may not be so large. The 
observed saturated ferromagnetism maintained under the 
gap formation seems to be related to such situations. 



G. Negative charge-transfer gap 

We here point out that the present material K^CrgOie 
is a doped negative charge-transfer-gap (CT-gap) type 
system^ in the Zaanen-Sawatzky- Allen phase diagram^ 
Because K^CrgOie is in the mixed-valent state, we need 
not invoke the negative CT-gap situation 18 for metalliza- 
tion of the system, as is realized in Cr02 where the metal- 
lization via self-doping mechanism operates^ However, 
in a hypothetical KoCrgOie system, where K ions are 
completely depleted and all the Cr ions are in the d 2 
state, this self-doping mechanism should be relevant for 
both metallization and double-exchange ferromagnetism 
unless the strong electron correlations leads the system 
to a Mott-insulating ground state with spin S — 1 lo- 
cal moments aligned antiferromagnetically or with local 
spin-singlet pairs of S = 1 spins. We here point out 
that, even in K^CrgOie, the negative CT-gap situation 
is realized between the 0(2) 2p and Cr 3d t2 g orbitals, 
as is evident in the calculated negative spin-polarization 
of the 0(2) ions. In fact, results of our GGA+U calcu- 
lation for K2Cr80i6 indicate that the gap in the DOS 
opens when U is large enough although the Fermi level 
is away from the gap due to the mixed-valent nature of 
the system. Thus, the deficiency of K ions, if realized 
experimentally, can offer an important opportunity for 
studying the negative CT-gap situation and its carrier 
number dependence. 
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IV. SUMMARY AND PROSPECTS 

In summary, we have made the first-principles elec- 
tronic structure calculations and predicted that a 
chromium oxide K^CrgOie of hollandite type should be 
a half- metallic ferromagnet. We have shown that the 
double-exchange mechanism is responsible for the ob- 
served saturated ferromagnetism. We have discussed the 
possible scenarios of the metal-insulator transition ob- 
served at low temperature, which include possibilities 
of the formation of the charge order and semimetal-to- 
semiconductor transition, but we have argued that the 
formation of the incommensurate, long- wavelength den- 
sity wave of spinless fermions caused by the Fermi-surface 
nesting may be the origin of the opening of the charge 
gap. We hope that these predictions will be checked by 
further experimental studies. 

Before closing this paper, let us discuss some prospects 
for future studies of this material, which may include the 
following: 

(i) This material offers an interesting opportunity for 
studying the effects of the majority-spin band gap on 
the electronic properties of half-metallic ferromagnets. 
In particular, the so-called non-quasiparticle state o 21 i 22 
and the effects of the opening of the gap on the non- 
quasiparticle states should be pursued by the spectro- 
scopic experiments such as photoemission and X-ray ab- 
sorption experiment j 23 i 24 i 25 i 26 

(ii) The transport properties of K^CrgOie are also 
anomalous, 7 as those of a "bad metal" ? 22^27^28^29^ which 
may also be interesting because not only the spin fluctu- 
ations but also the collective fluctuations of the orbital 
degrees of freedom between the degenerate d yz and d zx 
orbit als may play an important role. 

(iii) Introduction of the deficiency of K ions, if it could 
be made experimentally, would enable us to examine the 



effects of changing the doping rate on the gap forma- 
tion, double exchange ferromagnetism, and half metallic- 
ity. We also point out that only the ID-like Fermi sur- 
face should remain by ~10% hole doping if we assume 
the rigid-band shift of the Fermi level (see Fig. 5), sug- 
gesting the realization of the ID electron system, 
(iv) Theoretically, one should study the three-band Hub- 
bard model with the d xy , d yz , and d zx orbitals to see, in 
particular, the effects of orbital fluctuations on the elec- 
tronic states of the system. The study may be extended 
to the d-p model including the 2p orbitals of O ions if one 
wants to take into account the negative CT-gap situation 
of the present system. Techniques beyond the GGA+Z7 
method, such as the dynamical mean- field theory^ and 
variational cluster approach^ 2 - are required to examine 
the presence of non-quasiparticle states and their defor- 
mation due to the gap formation. 

Finally, we want to point out that, because the basic 
electronic states of the present system have many aspects 
in common with those of Cr02, it may be very useful to 
clarify similarities and differences in the electronic states 
between K^CrgOie and Cr02 in all the above respects. 
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